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两 两 NQD 序列 下 线性 指数 分 布 参数 的 经 验 Bayes 检验 * 


EOS XE 
(西北 工业 大 学 应 用 数学 系 ， 西 安 710072) 


摘 要 : 本 文 讨论 了 线性 指数 分 布 参数 的 经 验 Bayes 单 边 检验 问题 。 对 于 两 两 NQD 样本 序列 ， 利 用 核 函 
数 估计 方法 给 出 了 检验 函数 。 证 明了 该 检验 函数 是 渐 近 最 优 的 ， 并 在 适当 的 条 件 下 得 到 其 收敛 速 


度 。 
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1 引言 


经 验 Bayes (Empirical Bayes, EB) 方法 是 Robbins 于 上 世纪 五 六 十 年 代 提 出 的 。 对 于 EB 估 
计 和 检验 问题 已 有 了 相当 多 的 研究 。 考 虑 到 在 实际 问题 中 样本 大 都 存在 一 定 的 相依 关系 及 近年 
来 线性 指数 分 布 在 新 产品 ， 新 工艺 有 效 性 检验 等 方面 的 重要 意义 ， 本 文 在 两 两 NQD (Negative 
Quadrant Dependent, NQD) 序列 下 ， 讨 论 线性 指数 模型 参数 的 EB 检验 问题 ， 利 用 概率 密度 
核 估计 方法 给 出 参数 的 EB 检验 函数 ， 获 得 了 该 检验 函数 的 渐 近 最 优 性 及 收敛 速度 。 

首先 给 出 两 两 NQD 序列 的 定义 。 

定义 1 称 随 机 变量 (7.v.)X 和 YY 是 NQD 的 ， 若 对 任意 的 z, y € R， 有 


P(X <z,Y <y) < P(X <7)P(Y <y), 


Prv. EI (X4, :n = 1,2,…} 是 两 两 NQD IJ, GEHE Ej X; 与 Xi 是 NQD 的 。 
两 两 NQD 序列 是 由 Lemmanil 提出 的 ， 可 以 看 出 它 一 类 应 用 十 分 广泛 的 相依 序列 ， 在 可 靠 
性 理论 ， 多 元 分 析 等 领域 中 都 十 分 常见 ， 许 多 负 相 依 序列 都 是 在 它 基 础 上 生成 的 ， 如 著名 的 负 
相依 (NA) 序列 就 是 它 的 特殊 情况 之 一 ， 因 此 对 于 两 两 NQD 序列 的 讨论 就 显得 更 加 重要 。 在 这 
方面 已 有 了 一 些 研究 ， 如 Mutulal 引 获得 了 同 分 布 的 两 两 NQD 序列 的 Kolmogorov 型 强大 数 定 
律 。 
考虑 线性 指数 模型 ， 设 随机 变量 XX 的 条 件 密度 函数 为 


f(z10) = (uz + 0)exp[—07 — uz?/2], (1) 
其 中 60, 为 参数 且 j 已 知 ， 样 本 空间 Q = {r:r > 0}， 参 数 空间 
OQ -—40 : 0)dr = 1). 
{ >0 /relodz 1} 
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假设 参数 0 的 先 验 分 布 为 G(O)， 则 区 的 边缘 密度 及 其 导数 为 
OE / f(z10)dG(0), f'(z) = f (u — ja? — 2ubx — 0? )e- 05-97 Pag). 
e e 
考虑 单 边 检验 问题 ，Ho : 9 < bo Hi:0» bo， 其 中 b 是 给 定常 数 。 取 线性 损失 函数 
L(0,do) = max (c(0 -bo);0}， 工 (0,d) = max (c(69 — 0),0), 


这 里 c 为 正常 数 ，D = (do, di) 是 行动 空间 ，d; 表示 接受 H;, i = 0,1. 
取 随 机 化 判决 函数 5(z) = P{ 接 受 Ho |X =r} M (e) 的 风险 函数 为 


R(6(2),G) = 人 f [L(0, do)&(z) + L(0, d1)(1 — 9(2))] f (æ | 0)dzaG(0) 


= [sos + Co, 
这 里 
Cs - L L(0,d;dG(0), 
a(z) = js «(8 — bo) f (210)dG(0) = u(z) f(z) + v(x)f (x) + w(z)pc(2), 
其 中 


us) = -eur +0), vo) = —e wis) pel®) = [ oaG(0) 
eo 
进而 0 的 EB 检验 函数 为 ic = 1. 1a(z)<o)， 相 应 的 Bayes 风险 为 
Rlic,G) = int R(z),G) = f Bd (2) 
*4 G(0) 已 知 且 5(z) = jc(z) 时，Bayes 风 险 R(6G,G) 可 以 取 到 ， 但 通常 G(b) 未 知 ， 所 以 检 
验 函 数 c(z) 不 可 用 ， 从 而 引入 EB 方法 。 
2 两 两 NQD 序列 经 验 Bayes 检验 函数 的 构造 


设 (X1,01),… , (Xn,0n) 为 同 分 布 两 两 NQD 序 列 。 本 文 假定 f(z) € Cs,a， 这 里 Cs。 表 
示 呈 中 的 一 类 s 阶 导数 存在 、 连 续 且 绝对 值 不 超过 a 的 概率 密度 函数 族 ，s > 3 为 正 整数 。 

以 下 逐步 给 出 检验 函数 w(z) 估计 量 的 构造 。 

设 Klx) (i = 0,1) 是 区 间 (0,1) 上 有 界 可 测 函 数 ， 几 乎 处 处 可 微 ， 在 (0,1) 之 外 为 零 且 满足 


U l1, t=%, 
(L) — | wKi;(u)du- 
t! ; 
t 40 0, tzi, t=0,1. 


下 面 给 出 概率 密度 函数 f(x) 及 其 导数 f^ (a0) 的 核 函数 估计 为 


fala) = a. i), fne a: E) 
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其 中 {hn} ERUF], H hn — 0(n — oo). 
同时 ， 由 于 
十 co 
pe@) = fdr = E [Toon], 
H pa (x) 的 估计 为 . 
Pa(z) = LS aos 
i=l 
从 而 检验 函数 w(z) 的 估计 量 an(z) 为 an(z) = ult) fale) + u(r) fl) + w(x)pa(z). HN 


的 EB 检 验 函 数 定义 为 bo(z) = 1-165,()«o)« 9 Es RII X1, X2, ,Xn 的 联合 分 布 取 期 
望 ， 则 bu(z) 的 全 面 Bayes 风险 为 


Ru G) = E.| f al(z)5n(z)dz + Cc] = f o(z) E, [5 (z)]dz + Co. (3) 
a Q 
EX2 设 多 表示 参数 0 的 先 验 分 布 族 ， 若 对 任意 的 G(O) EF, di 
Jim R(,, G) = Ric,G)， 
则 称 EB 检 验 函 数 {in(z)} 为 渐 近 最 优 的 ; 车 R(6%,G) — R(ós,G) = O(n) (p > 0), W 
称 EB 检验 函数 具有 收敛 速度 O(n-?)。 
3 ”若干 引 理 


本 节 给 出 若干 引 理 。 这 里 令 M, M, 表示 不 同 常数 ， 即 使 同一 表达 式 中 也 可 取 不 同 值 。 

5| 10] RENZE X FLY 是 NQD Bf, #r, s 同时 为 非 降 ( 非 增 ) 实 函数 ， 则 r(X)，s(Y) 
仍 为 NQD 的 。 

由 定义 知 K(x) (i = 0, 1) 是 实数 域 尽 上 的 有 界 变 差 函数 ， 根 据 实 变 函 数理 论 ， 存 在 RR 上 的 
单调 有 界 变 差 函数 Ki (m), K2(x), i = 0,1， 使 得 Ki(z) = Ki(x) — K2(x), i=0,1。 

因为 {Xi :i = 1,2,… n] 为 两 两 NQD 序列， 由 引 理 1 知 


— Xir 
aR 





为 两 两 NQD 序列 。 进 一 步 有 


1 77 一 用 1 7 i : n 
G5. ( T hya $ (Qus i=1,2, m=0, {ou>} 5 {0} 


为 同 分 布 两 两 NQD 序列 ， 从 而 
1 n 1 n n n n 
f(s) = D0- 0, L= EG- ER, Pe ER, 
n j= D jz n jz nj n ja 
BI an(z) 的 各 分 量 可 以 分 解 为 两 两 NQD 序列 和 的 形式 。 
539 201 设 {Xi := 1,2……，,m} 为 两 两 NQD 序列 
jk 


EX? < +o, Ti(k)= M, QG- EX), j21 
i=j+1 
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则 
了 十 天 
E(T(k)) < V, EX?. 
?一 7 十 1 


引 理 3 向 4 R(ôc,G), R(ôn, G) 如 前 式 (2),(3) 定义 ， 则 

0< R(64,G) — R(ég,G) < c y) la(z)IP las (z) — a(z)| > la(z)l)dz 

QR 
384 设 {Xi:i=1,2,… n} ARA NQDJY JJ, BRAR (L) 满足 ， 则 
1) 3in—0, h,—0, nh}, > +o}, Ej|fa(z) — f(z)? — 0: 
2) 3*n0, hn — 0, nh > +t, E.|fz(x) — f'(z)? — 0; 
3) 34n08B. E,.pc(z)-pc(x)P? — 0. 
证 明 1) 由 Cr 不 等 式 
En|fn (7) 二 Fa? < MiE, [fa (z) e E, f. (x)]- + M3 [Es fa (2) 一 HOE 


一 方面 
IE. f.) - Fæ) = | | KOFE- thn) — Feat, 


因为 f(r)€ Css. Bi Taylor ÆRA 
fæ- tha) = f) + LE ane) + 


TEC thn) P Cay, E € (£ thua), 





利用 核 函数 性 质 可 知 
(Enh) - ft) < M| [ KEE rasa 





< Mhs. 
另 一 方面 ， 由 引 理 2 可知 
En [fn(x) — E, f(x J? < 2» {> Di (QG — EQG) y 


2 n 
< 2 y: Y ptt - B00) < 2p[oh] + EaR] 


E 

因此 结论 1) 成 立 且 l 
|E fs (2) - f(2)] < Mhay En fala) - Enfal)? < zz 

2) 类似 证 明 1)， 结 论 2) 成 立 且 
(Bf) - FO) s MA, E faz) - Esa < hr: 

3) 类似 证 明 1)， 结 论 3) 成 立 且 


M 


|Eno(2) - pao(z)| =0, Ea|pa(z) - Enele)? < 元， 
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4 ”经 验 Bayes 检验 函数 的 渐 近 最 优 性 及 收敛 速 度 


本 节 给 出 文章 主要 结果 ， 即 EB 检验 函数 的 渐 近 最 优 性 及 其 收敛 速度 。 
定理 1 ( 渐 近 最 优 性 ) 设 {X; :1= 1,2,… ,n} 为 两 两 NQD 序列 ， 核 函数 条 件 (L) 满足 ， 且 


f lelàG(6) < +00, 
e 
则 当 n — oo, hn — 0, nhá 一 co 时 ， 有 
dim R(ôn,G) = R(60, G). 
证 明 令 
D, (z) = la(z)lP(lan(z) — o(z)] > Jo(z)]). 
则 Dn(z) < lo(z)]. Hals) ZX, A 
f aljiri f IldG(0) + cgo < +00, 
Q e 
因此 由 控制 收敛 定理 及 引 理 3 得 
lim R(5,,G) — Róg,G) & c f iin Dalejdz: 
noo Q ^99 
由 Markov 不 等 式 及 a(z) 定义 有 
D«(z) < |u(z)|Es|fa(z) — f(z)) + |v(z) fi (x) — f'(x)] + |w(z)|EsIpe(x) -pc(z)|， 
根据 定理 条 件 及 引 理 4 得 ， 对 任意 的 x & 9， 有 lim Dn(z) = 0， 即 
lim R(5,G) = Ric,G) 


定理 2 (收敛 速度 ) 设 {X;:i= 1,2,… ,n} 为 两 两 NQD FJ], RAAR (L) WE. HX 
任意 的 入 e (0, 1) 
Í | dz « oo, f aal ^ dz « oo, 
Q Q 
则 当 h = n- En, A 
R(6,,G) — R(ôc, G) = o(n- 3658). 
证 明 由 Markov 不 等 式 及 引 理 3 可 知 
R(&,G) — RG6G,G) € MiEs|fa(z) — f(2)) + MoEnlfs(z) — f(z) 
*M3Es [Pa (2) ~ pe(z)|, 
由 条 件 0 < 入 <1<2 及 引 理 4 可 得 
En|fn(z) - F(E) < [Enslfa(z) - Fœ) 


[Mi [Es fa (2) — f(x)? + M2Enlfn(z) — En fao ^ 


I^ 


I^ 


? 


1- nh2**? A/2 


m| nh2 
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类 似 可 得 
2s 十 2 入 /2 
Enlia) - GI s [RM 


A/2 
nhá : 


1 
,Bnlpc(z) - pe) < m| =] 
k. “hrs mn 时， 有 


R(ön,G) — R(ôc,G) < a et E P ag 


nh2 nhá 


A(2s—5 
< Mn 3089), 


n 


A(2s—5) 
R(ôn, G) — R(ôc, G) = O( n-as ) 


Y 
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Empirical Bayes Test for Parameters of the Linear Exponential 
Distribution with NQD Random Samples 
WANG Liang, SHI Yi-min 


(Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072) 
Abstract: In this paper, the one-sided test problem for parameters of the linear exponential families 
was investigated through the empirical Bayes (EB) approach. The EB test function was constructed 
by using the kernel-type estimation for Negative Quadrant Dependent (NQD) samples, and the con- 
vergence rates for the EB test function was obtained under suitable conditions. The results show that 
the EB test function is asymptotically optimal. 

Keywords: NQD samples; kernel estimation; empirical Bayes test; asymptotic optimality; convergence 
rates 


Received: 04 Nov 2008. Accepted: 24 Dec 2009. 
Foundation item: The National Natural Science Foundation of China (70471057); the Natural Science 
Foundation of the Education Department of Shaanxi Province (03JK065). 





